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Nucleation of Vortex State in Ru-inclusion in Eutectic Ruthenium Oxide Sr 2 Ru0 4 -Ru 
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Eutectic samples of S^RuCU-Ru possess many /^m-sized Ru-metal inclusions. Based on a 
Ginzburg-Landau formulation we analyze the interplay of the chiral p-wave state of Sr2RuC>4 
(T c = 1.5-K") and the s-wave state of Ru metal (T c ,Ru = 0.5K) for an inclusion of cylindrical 
geometry. As a consequence of the mismatch of the order parameter phase, the occurrence of 
a spontaneous flux distribution appears for T < T c ,r u which evolves into a "Josephson vortex" 
on the Sr2Ru04-Ru interface. At a sufficiently low temperature a depinning transition can 
occur whereby the vortex moves to the center of the cylinder. 
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Sr2Ru04 is considered a textbook case of an uncon- 
ventional superconductor. Over many years evidence has 
accumulated that its pairing symmetry is of chiral p- 
wave type, i.e. spin-triplet pairing with broken time re- 
versal symmetry, generally denoted by the vector d = 
z(k x ± iky) (p x ± ip y )}~ 3 The critical temperature is 
strongly sensitive to impurities and reaches T c — 1.5K 
in the cleanest samples. It came as surprise when Maeno 
and coworkers discovered that in eutectic samples of 
Sr 2 RuC>4 with /zm-sized Ru-metal inclusions the onset 
of inhomogeneous superconductivity could be found at a 
temperature roughly twice that of the bulk critical tem- 
perature. 4 Consequently, this superconducting phase was 
dubbed "3-Kelvin phase" (3K-phase). It was early spec- 
ulated that the 3K-phase had filamentary nature and 
originates from nucleation of superconductivity at the in- 
terface between Ru-metal and Sr2Ru04. 5 Ru-metal itself 
is a conventional superconductor with at T c r u = 0.5K. 
It has not so far been possible to identify the micro- 
scopic mechanism which causes the increased transition 
temperature of the 3K-phase at the interface. 

Unlike in conventional superconductors where inhomo- 
geneous nucleation of superconductivity leads through 
simple percolation eventually to a bulk superconducting 
phase, it was suggested that the unconventional nature of 
Cooper pairing yields a more complex evolution from the 
3K-phase to bulk superconductivity in Sr2Ru04 as tem- 
perature is lowered. 5 ' 6 Indeed there is even a symmetry 
breaking transition on the way to bulk superconductiv- 
ity, since the nucleation of a p-wave pairing state on the 
Sr2Ru04-Ru interface leads to a time reversal conserv- 
ing phase incompatible with the bulk superconducting 
phase. This second transition can actually be identified 
in experiment, e.g. in anomalies in the critical current 7 
and in quasiparticle tunneling 8-10 as was shown in Ref. 6 

So far the situation for these eutectic sample at low 
temperature, where also the Ru-inclusions become su- 
perconducting on their own, has not drawn much atten- 
tion. In this study we would like to address one particular 
aspect of this situation which is connected with the fact 
that actually the two superconducting phases, s-wave for 
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Fig. 1. Flux distribution and vortex state: (a) immediately below 
T C: r u a sinusoidal flux pattern appears with two nodes diametral 
on cylindrical interface; (b) at lower temperature gradually a 
more concentrated flux line emerges on the interface; (a) and (b) 
correspond to a state of ip without phase winding; (c) below a 
critical temperature the flux line depins from the interface and 
forms a vortex at the center of the cylinder, a state with phase 
winding of rjj. 
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Ru and chiral p-wave for Sr 2 Ru04, meeting at interface 
are not phase compatible, a feature we will explain be- 
low. We would like to demonstrate here that this could 
lead to the creation of intriguing magnetic flux distribu- 
tions which may undergo a transition between two dis- 
tinct states as temperature is lowered. 

We illustrate this idea using the convenient geome- 
try of a cylindrically shaped Ru-metal inclusion in an 
infinitely large matrix of Sr 2 Ru04, which is most suit- 
able for the chiral p-wave phase. The axis of the cylinder 
lies along the z-axis and the radius R shall be of order 
of several (im, a dimension we assume to be large com- 
pared coherence lengths of both superconductors. When 
Ru becomes superconducting the p-wave superconduct- 
ing state of Sr2RuC>4 is already solidly established and 
will be considered as rigid here. Looking now at the su- 
perconducting phase of Ru we may model it by the fol- 
lowing Ginzburg-Landau (GL) free energy. We choose 
here cylindrical coordinates (r,6,z): 
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where F denotes the free energy per unit length along 
the z-axis with real coefficients a = a' (T — T c ,r u ), b and 
K and the covariant derivative D = V — i2eA/hc (ne- 
glecting z-dependence). The first part is the standard 
GL free energy for the order parameter ip of the conven- 
tional superconducting phase within the Ru-inclusion. 
The second term is for our purpose most essential as 
it represents the coupling between the order parame- 
ters for the two materials with the coupling constant 
K. At this level it describes the lowest-order Joseph- 
son coupling. 11, 12 For the chiral p-wave state we define 
ff = (r] x ,f]y,0) = rj p (l,±i,0) with a complex constant 
amplitude r] p and n as the normal unit vector of the in- 
terface. Note that for this Josephson coupling spin-orbit 
coupling is important as it provides the connection be- 
tween the spin-singlet and spin-triplet configurations of 
the Ru s-wave and the Sr 2 Ru04 p-wave phase, respec- 
tively. The structure of the coupling incorporates the se- 
lection rule of conserved total angular momentum J„ of 
the Cooper pairs along the interface normal (J ■ n = J n 
conserved). 11, 12 

The interface term imposes boundary conditions for ip 
at r = R, 



Kn ■ Dip — Kz ■ (ft x if) 

r=R r=R 



(2) 



Taking the chiral p-wave component as rigid and n — 
(cos 9, sin 6, 0) we may write for the term on the right 
hand side, 



z ■ (ft x ff) = —irj p e %e = \r] p \e 



(6-a) 



(3) 



choosing among the two chiral phases, the p x — ip y state 
(with a gauge transformation we may choose a — which 
is equivalent to the rotation of the coordinate frame 
around the z-axis by the angle —a). The resulting bound- 
ary condition implies the competition between two low 



energy configurations for ip. One option is a state with- 
out angular dependence of ip which ignores the coupling 
to the phase imposed by the chiral p-wave state through 
the interface. This state is stabilized at the expense of 
interface energy. The other is to pick up the phase wind- 
ing of the chiral p-wave state with the disadvantage of 
having to introduce an energetically costly singular line 
(vortex) which for symmetry reason would lie on the axis 
of the cylinder (r = 0). If the coupling at the interface is 
weak, it is the former phase which nucleates at T = T c ,r u 
and is stable for some temperature range. At lower tem- 
perature the other phase becomes competitive and may 
eventually win. 

In order to obtain a qualitative understanding of the 
evolution of these states we focus now on the interface 
assuming the order parameter ip to be finite and rigid in 
its modulus, such that we may restrict to discussion of a 
variable order parameter phase, 



iP(r = R,6) = \iP{R)\e im 



(4) 



at r = R. Then we can derive in a standard way the 
effective free energy functional for the phase <p(9) at the 
interface which has the form 

- A 2 cos(0 - 6) (5) 
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with the energy scale 
fo = 



2iTy/d e ffR 
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and A 2 / = 2^|^||r/j,| where d e ff = X s + X p + d is the ef- 
fective width of the interface for magnetic fields including 
the two London penetration depths, A SjP and the width 
of the interface region, d. 13. The temperature dependence 
of A 2 is governed by the order parameter \ip(T)\ and the 
London penetration depth A s (T) . We consider the other 
quantities such as \r) p \ and \ p as basically independent of 
temperature. Then the evolution as temperature is low- 
ered corresponds to a monotonic increase of A 2 from zero 
at T = T c ^r u . We introduce here also the magnetic flux 
quantum $ = hc/2e of a superconductor. 

Let us now discuss the behavior of the phase 4>{9) at 
the interface in our simplified model. The variation of 
Eq.(5) with respect to <p(9) leads to 



d 2 <P 
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Note that the derivative of (p(9) with respect to 9 corre- 
sponds to the local magnetic field in z-direction, 



B(9) = 



-d(p{9) 
1 09 



(8) 



like in standard Josephson junctions. 13 We turn now the 
solution of the variational equation. Within the given 
boundary condition that (p(9) = <p(9 + 2ir) without any 
phase winding for the order parameter ip. We consider 
approximative solutions for two limiting cases, 

A 2 sin 9 for A < 1 , 

m « { (9) 

9-4 arctanfe^-^] for A > 1 . 
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Fig. 2. Free energy balance for the interface state: for small values 
of A ("high" temperature) the state of ip without phase winding 
(flux distribution restricted to the interface) is realized; for large 
enough A the vortex state with phase winding of ip is stabilized. 
The transition can be viewed as a spontaneous flux line depinning 
transition. 



which yields to the magnetic field distribution 

A 2 cos 9 for A < 1 
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for A > 1 



(10) 

These solutions show a magnetic flux distribution lo- 
cated around the interface. The net flux vanishes due to 
flux quantization in the superconductor, since there is no 
phase winding of the order parameter of S^RuO.^ With 
decreasing temperature the flux distribution changes 
from a soft sinusoidal form to a more and more uneven 
structure where one magnetic field direction becomes 
spatially more concentrated evolving into a Josephson 
vortex, and the magnetic flux in opposite direction is 
more spread (Fig.l a,b). Also the local magnetic fields 
grow with shrinking temperature. Note that in our model 
geometry there is a rotational degeneracy for this flux 
distribution. 

While the state with this topology is stable for some 
range of temperature below T c ^r u , the alternative state 
with an order parameter i/j(r, 9) = \ifj(r)\e possessing a 
single vortex at the center of the cylinder becomes en- 
ergetically more competitive at lower temperature. The 
energy of this vortex state can be estimated as 



F v = 2nf 
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where the first term is obtained from Eq.(5) by setting 
<j)(9) = 9. The second term is the line energy per unit 
length of the vortex which we approximate by 
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with k = A/£ and e c the core energy which is a fraction 
of the first term constituting the magnetic contribution 



to the line energy. 

In Fig. 2 we show the energy F = F int - 2tt/o(1/2 - A 2 ) 
obtained by numerical solutions of the Sine-Gordon equa- 
tion in Eq.(7). For A > 1 the energy follows an essential 
linear behavior: F/fo ~ —7.2 + 8.2A. The transition be- 
tween the two states is determined by F ~ e„, i.e. the 
condition when the line energy of the vortex at the center 
of the Ru-inclusion and the energy to carry the flux dis- 
tribution on the interface are equal. Neglecting the core 
energy e c , taking InK ~ 1 and d e ff ~ 2A, we obtain the 
rough criterion, 



F(A) 
fo 
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Assuming a low-temperature value of A ~ 0.2/im and 
R ~ 10/im we estimate from Fig. 2 that the criterion is 
satisfied for A w 3.2. It is obvious that with increasing 
R the critical value of A for the transition grows. Con- 
sequently, for large Ru-inclusions the transition between 
the two states is less likely to occur than for smaller in- 
clusions. The state after the transition corresponds to 
a single vortex at the center of the Ru-inclusion with a 
compensating magnetic flux on the interface to yield an 
overall vanishing flux (see Fig.lc). 

In summary our study shows that the coupling of 
the intrinsic superconducting order parameter of a Ru- 
inclusion in the eutectic ruthenate samples is subject to a 
frustrating coupling with the p-wave order parameter of 
the surrounding Sr2Ru04 inducing states with a sponta- 
neous magnetic flux distribution. There are two compet- 
ing states which both induce a spontaneous flux pattern: 
the interface vortex state and the Ru-center vortex state, 
which are distinguished by the topology of the order pa- 
rameter ip (phase winding) of the Ru-superconductor. 
We show here that there is a clear sequence, how these 
states would appear as temperature is lowered. First the 
interface vortex phase appears below T Cj # u and only at 
lower temperature a transition to the center vortex state 
occurs. This transition can be viewed as a depinning 
transition of the interface vortex which moves to the 
center of the cylinder and leaves behind a compensat- 
ing counter flux uniformly spread over the interface. This 
suggests also that the additional pinning effects, through 
interface defects or special conditions due to varying in- 
terface curvature, would play an essential role too in 
defining the transition point. Such issues are certainly 
important for inclusions of rather irregular shape. Even 
under these condition transitions between different stable 
flux configurations are possible, although the situation is 
considerably more complex and beyond the scope of our 
study. We believe that these effects should be accessible 
to studies using scanning Hall or SQUID probes. 
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